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Abstract. We propose a new hypothesis for the origin of 
the major part of non-solar hadronic cosmic rays (CRs) 
at all energies: highly relativistic, narrowly coUimated jets 
from the birth or collapse of neutron stars (NSs) in our 
Galaxy accelerate ambient disk and halo matter to CR en- 
ergies and disperse it in hot spots which they fornr when 
they stop in the Galactic halo. Such events are seen as 
cosmological gamma-ray bursts (GRBs) in other galaxies 
when their beamed radiation happens to point towards 
Earth. This source of CRs is located in the Galactic halo. 
It therefore explains the absence of the Greisen-Zatsepin- 
Kuz'min cutoff in the spectrum of the ultra-high energy 
CRs. The position in energy of the "ankle" in the CR 
energy spectrum is shown to arise in a natural way. More- 
over, an origin of lower energy CRs in the Galactic halo 
naturally accounts for the high degree of isotropy of CRs 
around 100 TeV from airshower observations, and the 
small galactocentric gradient of low-energy CRs derived 
from ganrma-ray observations. 
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1. Introduction 

The origin of non-solar cosmic rays, discovered in 1912, 
is a major unsolved puzzle in physics. Recent observa- 
tions do not confirm some expectations in the currently 
accepted framework for hadronic CR origin, where cosmic- 
ray nuclei with energies below 3 x 10^^ eV (the "knee") 
are accelerated in Galactic supernova remnants (SNRs) 
(Ginzburg 1957), and those above 3 x 10^^ eV (the "an- 
kle") (UHE CRs), for which a disk origin is unlikely 
due to their isotropy, in sources far beyond our Galaxy 
(Burbidge 1962). Some observational anomalies are: 
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— The "GZK cutoff" in the CR intensity at energies 
above ^ 10^" eV, due to interactions of the extragalac- 
tic CRs with the cosmic microwave radiation, as pre- 
dicted by Greisen (1966) and by Zatsepin & Kuz'min 
(1966) seems to be absent (Takeda et al. 1998). At the 
same time the directions of the events with the highest 
energy do not line up with nearby potential extragalac- 
tic CR sources. It is difficult to understand these two 
observations within any model of an extragalactical 
origin of UHE CRs (HiUas 1998). 

— No very-high energy 7-rays, due to proton in- 
teractions, are observed from some SNRs, for 
which a measurable flux is expected in simple 
models (Aharonian et al. 1994) within the conven- 
tional framework (Prosch et al. 1996; HeB et al. 1997; 
Buckley et al. 1998). 

— Detailed models (e.g., (Strong & Moskalenko 1998)) 
yield a significantly larger galactocentric gradient in 
the sky distribution of high-energy (>100 MeV) 7- 
rays from interaction of CRs from SNRs in the Galac- 
tic interstellar medium than the one observed by the 
EGRET detector on the Compton Gamma-Ray Ob- 
servatory (Hunter et al. 1997; Strong & Mattox 1996; 
Erlykinet al. 1996). 

— Diffusive propagation of CRs from sources sharply 
concentrated towards the Galactic disk and centre 
(such as SNRs) yields anisotropics in the distribu- 
tion of charged CRs at an energy of about 100 TeV 
in excess of the observed value (Aglietta et al. 1995) 
by more than an order of magnitude. This is consis- 
tently found both in simple order-of-magnitude esti- 
mates (Gaisser et al. 1995) and sophisticated calcula- 
tions of CR propagation (e.g., Ptuskin et al. 1997). 

The last three observations suggest that, perhaps, SNRs 
are not the main accelerators of Galactic hadronic CRs at 
energies below the knee. 

Maybe these anomalies can be resolved in more sophisti- 
cated models within the conventional framework; e.g. by 
modifying standard ideas about intergalactic CR propa- 



A.Dar & R.Plaga: An alternative source of cosmic rays at all energies 



gation for the first (Sigl ct al. 1998) and interstellar CR 
propagation for the third and fourth (Erlykin et al. 1997) 
anomaly, respectively. Spectral cutoffs at energies much 
below the knee could be the reason for the absence of TeV 
7-rays from some fraction of SNRs (Baring et al. 1997). 
Here we explore the possibility that the anomalies are 
first indications that the current framework for the ori- 
gin of hadronic CRs is wrong, and propose an alternative 
origin. 

Gamma-ray bursts (GRBs) are short and intense bursts 
of 7-rays from distant galaxies that are detected by 
space satellites at a rate of -Rgrb ^ 10^ per year 
(Fishman & Hartmann 1997). Their origin is unknown, 
but various observations suggest an association with the 
formation of black holes or the birth or collapse of neu- 
tron stars (NSs) due to mass accretion or phase transition, 
in close binaries or alone. The strongest indication for a 
connection with the birth of NSs are the coincidences of 
supernovae with GRBs, e.g. SN 1998bw and GRB 980425 
(Galama et al. 1998) and of SN 1999E and GRB 980919 
(Thorsett & Hogg 1999; Kulkarni & Frail 1999). A con- 
nection between events in the life of NSs and GRBs at 
cosmological distances was first suggested by Paczynski 
(1986), Goodman et al. (1987) and Dar et al. (1992). 
The rest of the introduction discusses the assumptions we 
make about GRBs in our scenario, beaming in GRBs and 
previous work on CR acceleration in GRBs. Sect. 2 turns 
to a more detailed quantitative discussion of our picture 
of GRBs. Sect. 3 discusses the propagation of the ejecta of 
GRBs in the Galaxy and Sect. 4 treats energy spectrum 
of CRs which is expected in our scenario. Sect. 5 and 6 
describe how our scenario resolves the above-mentioned 
anomalies of UHE CRs and lower energy CRs respectively. 
Sect. 7 concludes the paper. 

1.1. Basic assumptions for our scenario 

We make two basic assumptions for our scenario. 

— GRBs are not rare events, but common phenomena in 
the life of neutron stars; i.e., in our Galaxy they occur 
with a frequency comparable to the NS birth rate. We 
argue in the next subsection that a necessary implica- 
tion from this and the observed rate of GRBs i?GRB 
- that the emission from GRBs is strongly beamed - 
seems likely in view of recent GRB observations. 

— Jetted ultra-relativistic ejecta ( "plasmoids" ) emit the 
beamed radiation in GRBs and carry a considerable 
fraction of the total binding energy of a NS. The ejecta 
remain confined to a small size until slowed down - via 
interaction with an ambient medium - to mildly rel- 
ativistic speeds. Only then do they dissipate most of 
their kinetic energy - partly by particle acceleration - 
and release most of the previously accelerated parti- 
cles. We give arguments in favour of this assumption - 
based on a special interpretation of NS kicks and analo- 



gies with similar astrophysical systems - and describe 
it in more detail in Sect. 2. 

We return to the question whether these assumptions are 
overly speculative in the conclusion. 

1.2. Beaming in GRBs 

If the true rate of GRBs is similar to the birth rate 
of NSs, Rns — 0.02 ycar^^ in galaxies like our own 
(van den Bcrgh & Tamman 1991), the radiation emitted 
in GRBs must be narrowly beamed into a solid angle 
Art ~ TT X 10~^ in order to explain their rate inferred 
from i?GRB (Wijers et al. 1997) in such galaxies: 



ii'GRB(MW) ~ 2i?Ns(Al^/47r) ^ 10"* year" 



(1) 



Indeed GRBs could be produced by highly relativis- 
tic, collimated jets with bulk motion Lorentz fac- 
tors F = 1/^/1 -v^/c^ - 10^ (Shaviv & Dar 1995; 
Blackman et al. 1996; Meszaros & Rees 1997; 

Chiang & Dermer 1997; Dar 1998a). 

In this connection - and always in this paper - the term 
"jet" is meant to imply only collimated ejecta into a 
small opening angle and not an emission stationary in 
time. The very large redshifts recently measured in some 
GRBs and their host galaxies - e.g. z=3.42 for GRB 
971214 (Kulkarni ct al. 1998) - seem firmly to require jet- 
ted emission for energetic reasons. For instance, the as- 
sumption of isotropic emission implies that > 4 x 10^^ 
erg - which is hardly obtainable in the birth or collapse 
of any stellar object - was emitted by GRB 990123 as 
7-rays (Djorgovski et al 1999; Kippen et al. 1999 ), while 
only ~ 10^* erg was emitted if the emission was narrowly 
beamed into a solid angle Ail ^ n x 10^^. Hjorth et 
al.(1999) found a complete absence of optical polarisation 
in GRB 990123 and mention a relativistic jet which is 
seen at a small viewing angle as one possible explanation 
for their somewhat surprising observation. Emission from 
"narrow" jets with bulk motion Lorentz factors F ^ 10'^ is 
indeed beamed into a solid angle Ail ^ tt/F^ ~ tt x 10^® 
consistent ^ with Eq.(l). Here "narrow" means that the 
angle of the emitted ejecta 9 subtended from the burst site 
is smaller than 1/F. There is some evidence that under the 
assumption of jetted ejection this condition is fulfilled in 
GRBs (though the following cannot be considered as proof 
for jetted ejection): Schaefer & Walker (1999) argued from 
timing data of the B ATSE detector that for the burst GRB 
920229, < 1'; this fulfils the "narrow" -jet condition for 
F < 3000. The amount of beaming and physical structure 
of the GRB jet assumed in our paper is very similar to 
the one in a model of Chiang & Dermer (1997) for GRB 
afterglows. 



^ Such Lorentz factor values are required from GRB obser- 
vations independent of whether beaming takes place (see, e.g.. 
Baring & Harding 1997). 
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If such strong beaming is typical - and this we assume 
- we observe only a very small fraction of a large rate of 
the events that can produce GRBs. We call these events 
"Galactic" (if they occur in our Milky Way (MW) galaxy) 
and "Cosmological" (if they occur in distant galaxies) 
gamma-ray bursts" (GGRBs and CGRBs), respectively. 

1.3. CR acceleration in GRB ejecta 

Gamma-ray bursters have been proposed as efficient ac- 
celerators of ambient material to CR energies before, 
but GGRBs with isotropic CR emission (Dar et al. 1992; 
Milgrom & Usov 1996) cannot supply the bulk of the 
Galactic cosmic rays for energetic reasons. Nor can 
the idea of UHE CRs from CGRBs (Vietri 1995; 
Waxman 1995; Bottcher & Dermer 1998) explain the ab- 
sence of the GZK cutoff (Dar 1999). In this paper we show 
that the effects of strong beaming make GGRBs plausible 
sources of non-solar hadronic cosmic rays at all energies 
(from ~ 100 MeV to 3 x 10^" cV) observed near Earth. 



Fig. 1. A highly schematic sketch of our scenario. The birth or 
collapse of NSs in the disk of our Galaxy leads to an ejection 
of two opposite jets that produce "hot spots" when they stop 
in an extended Galactic halo. 



2. Quantitative properties of the plasmoids ejected 
in GRBs 

The ejection of highly relativistic jets from accreting 
or collapsing compact objects, and especially the physi- 
cal mechanism for their observed coUimation "remains a 
major unsolved problem in astrophysics (Longair 1997)" . 
Therefore, instead of relying on theoretical models, we 
have tried to estimate the properties of those jets that 
may produce the GRBs and cosmic rays from data about 
NSs and - by way of analogy - from observations of similar 
jets. 

NSs are observed to possess a large mean velocity of 
V « 450 ± 90 km s~^ (Lyne & Lorimer 1994). Burrows & 



Hayes (1996) argue that an anisotropic stellar collapse - 
rather than proper motions due to orbital speeds in bina- 
ries - is the mechanism responsible for the necessary "kick" 
and and that matter - rather than neutrinos - carries most 
of the momentum. Janka & Raffclt (1998) find that works 
in which deformed neutrino spheres impart large kicks to 
NSs are incorrect and that it is in general difficult to im- 
part the necessary momentum via neutrinos. In this light 
- and making the above assumptions - it seems plausible 
that the ejection of relativistic jets in the birth or collapse 
of NSs is responsible for their observed large mean veloc- 
ity. Momentum conservation then implies that the kinetic 
energy of the jets is 



E 



K 



cP ^ vMnsc ^ Ax W erg 



(2) 




where we used the typical observed NS mass, Mms = 
1.4Mq. If two antiparallel jets are ejected, then the above 
estimate becomes a lower limit for the jet kinetic energy, 
that may be of the order oi Ey^ ~ 10^^ erg. Thus, already a 
modest fraction / ~ 0.01 of the total energy injected into 
the MW in jets from GGRBs, at a rate R ~ 0.02 year"!, 
similar to the NS birth rate, can supply the estimated 
Galactic CR luminosity Lmw[CR] ~ 1-5 x 10"*^ erg (e.g., 
Drury et al. 1989). 

Highly relativistic jets are emitted by all astrophysical sys- 
tems where mass is accreted at a high rate from a disk onto 
a central black hole (BH). They are observed in Galac- 
tic super luminal sources, such as the micro-quasar GRS 
1915-f 105 (Mirabel & Rodriguez 1994), where mass is ac- 
creted onto a stellar BH, and in many active galactic nuclei 
(AGN), where mass is accreted onto a super- massive BH 
(Begelman et al. 1984). Mildly relativistic jets from mass 
accretion were seen both in AGN and in star binaries con- 
taining NSs such as SS433 (Margon 1984) and Sco X-1 
(Fomalont et al. 1983). 

High-resolution radio observations of the micro quasar 
GRS 1915-1-105 resolved the narrowly coUimated rel- 
ativistic jets into clouds of plasma (plasmoids) that 
are emitted in injection episodes which are correlated 
with sudden removal of the accretion disk material 
(Mirabel & Rodriguez 1994). In GRS 1915+105, during 
the first five days after ejection, these plasmoids appear 
to have expanded with the speed of sound in a relativistic 
gas, Cs ~ c/v3 in the plasmoid rest frame, corresponding 
to vertical expansion speed of ^ c/Fvo in the lab frame), 
until they reach a radius Rp ^ 2 x 10^^ cm (^ 0.001 pc) 
(Rodriguez & Mirabel 1998; Mirabel & Rodriguez 1994). 
Their expansion seems to be slowed afterwards, proba- 
bly due to some confinement mechanism 
(Mirabel & Rodriguez 1994; Mirabel & Rodriguez 1999). 
Likewise, jets from quasars and radio galaxies can also 
often be resolved into the ejection of distinct clouds at 
the pc scale (e.g., Giovannini et al. 1999) and might be 
governed by mechanisms similar to the jets emitted by 
micro quasars. These jets retain a constant radius af- 
ter initial expansion (see, e.g., (Swain et al. 1998)) and 
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some confinement mechanism is clearly at work. The re- 
peated ejection of plasmoids (often called "plasnions" in 
the older literature) is a well known model for cxtragalac- 
tic jets (Begelman et al. 1984; de Young & Axford 1967; 
Christiansen 1973). 

The confinement mechanism which collimates plasmoids 
is not understood; two possibilities are magnetic con- 
finement (Begelman et al. 1984) and inertial confine- 
ment via the ram pressure of the ambient medium 
(de Young & Axford 1967). The plasmoids slow down by 
sweeping up the ionised interstellar medium in front of 
them. 

It is one of our assumptions that the plasmoids from the 
birth/collapse of NSs evolve in a similar way to the ones 
apparently ejected in galactic superluminal sources and 
AGNs (Dar 1998a). The plausibility of this idea was re- 
cently endorsed by Mirabel & Rodriguez (1999) in their 
review about relativistic jets the Galaxy; they write: "the 
study of the less extreme collimated outflows in our own 
Galaxy may provide clues for a better understanding of 
the super-relativistic jets associated with ... GRBs". The 
plasmoids in GRS1915+105 could only be followed up to 
distances of about 0.03 pc from the central source. A non- 
thermal jet with a radius of about 0.1 pc at a distance of 
about 60 pc from the central source GRS1915-I-105 and 
pointing back to it was observed recently in the radio 
range (Rodriguez & Mirabel 1998a). No connection be- 
tween the central source and this jet was observed. A 
possible interpretation is that this jet is a collimated plas- 
moid that was ejected by the central source. These various 
observations of jets motivate the first part of our second 
assumption in the Introduction that: 

— after an initial expansion - and until their final dissipa- 
tion phase - the plasmoids retain a radius of about 0.01 
pc, collimated by the as yet ill-understood mechanism 
apparently at work in the plasmoids and jets discussed 
above. 

Plasmoids ejected in the birth/collapse of NSs - as opposed 
to the ones in micro quasars and AGNs - are "isolated", 
i.e. their central source is active on a time scale smaller 
then their propagation time. Our assumption is that this 
fact docs not invalidate the qualitative similarity between 
the mentioned three object classes. 

In extragalactic double-lobe radio galaxies like Cyg-A 
(Begelman et al. 1984) it is observed that apparently bulk 
kinetic energy of the jet is dissipated mainly in a "hot 
spot" . The nonthermal jet at a large distance from 
GRS1915-I-105, mentioned above, is an indication that 
"isolated" plasmoids might show a similar behaviour. A 
rising internal kinetic pressure - due to the sweep up of rel- 
ativistic particles - and dissipation of their internal mag- 
netic energy might lead finally to their expansion and stop- 
ping in "hot spots" , where a considerable part of their 
initial kinetic energy is released in highly relativistic par- 
ticles. Another factor leading to a final expansion could 



be a decrease in inertial confinement due to the a falling 
ram pressure in decelerating plasmoids, but the detailed 
physics of hot spots produced by "isolated" plasmoids re- 
mains unclear. This set of observations motivates the sec- 
ond part of the second assumption in the Introduction: 

— the plasmoids dissipate most of their energy in turbu- 
lent motion of inter-halo matter and CR production 
only when they have been slowed down - via sweep- 
ing up ambient matter - to mildly relativistic energies. 
Most of previously accelerated particles are also re- 
leased in the phase. 

In the case of a GGRB there is no stationary fiux of parti- 
cles whose energy is dissipated at a relatively sharply de- 
fined "working surface" , as in the case of the jets of active 
galaxies (Begelman et al. 1984). Under the above assump- 
tions most of the dissipation and particle acceleration will 
take place in a region which is small compared with the to- 
tal distance travelled by the plasmoid and which we might 
still call "hot spot" . 

The possibility that relativistic jets responsible for 
GRBs are sometimes formed in supernova (SN) explo- 
sions has been discussed by various authors (Gen 1998; 
Wang & Wheeler, 1998; Eichler & Levinson 1998). It is 
conceivable that a considerable fraction of SN explo- 
sions leads to GRBs with the properties required in our 
scenario. The jet like features observed in many SNRs 
(Gaensler et al. 1998) may even hint towards this possi- 
bility (though they certainly do not prove it). However, 
this idea is not mandatory for our scenario; other events 
in the life of NSs might be responsible for most GRBs. 
Because the accretion rates and magnetic fields in- 
volved in events leading to GRBs are larger than in all 
other systems discussed above, an initial Lorentz fac- 
tor of r ~ 10^, as required by GRBs observations, 
does not appear excessive. The initial equipartition mag- 
netic fields in such plasmoids can be estimated to be 
B > 100 G. Extragalactic jets, which might be qual- 
itatively similar objects, have been identified as pro- 
lific CR accelerators (Burbidge 1962), in particular their 
hot spots (Biermann & Rachen 1993). Highly magnetised 
relativistic plasmoids are both efficient CR accelerators 
(Chiang & Dermer 1997; Dar 1998b) (through Fermi ac- 
celeration) and strong emitters of beamed 7-rays through 
synchrotron emission, inverse Compton scattering and res- 
onance scattering of interstellar light. When they point in 
(or precess into) our direction from external galaxies, they 
produce the observed GRBs and their afterglows. 

3. Propagation of the ejected plasmoids in the 
Galaxy 

For our discussion, we adopt an extremely simplified pic- 
ture of our Galaxy, namely, a Galactic thin disk sur- 
rounded by a spherical halo with an approximate radius 
of i?h ~ 50 kpc ( Vh ~ 1.5 X 10™ cm'^) and a mean gas 
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density of rih ~ 10~^ cm^'^. The highly relativistic plas- 
moids which are emitted with Ek ^ 10^^ erg perpendic- 
ular to the Galactic disk are assumed to be decelerated 
mainly by sweeping up ambient matter. Whether this is 
an accurate description depends on whether the swept-up 
material "sticks" to the plasmoid and is incorporated into 
it, rather than e.g. forming a hot expanding wake behind 
the plasmoid. A full dynamic treatment of the propaga- 
tion of ultrarelativistic plasmoids is required, a problem 
that is however beyond the scope of this paper. The plas- 
moids then stop in the Galactic halo, when the rest mass 
energy of the swept-up ambient material is equal to their 
initial kinetic energy: The column density of gas required 
to stop jets with a radius smaller than Rp ~ 0.01 pc is 



N, > 2.2[{EJ10'" CTg)/{Rp/O.Ol pc)^)] x lO" 



(3) 



The mean column density of gas (mainly molecular, 
atomic and ionised hydrogen) perpendicular to the Galac- 
tic disk, is Nc < 10^^ cm^^ (Bcrezinskh et al. 1990). From 
Eq.(3) one reads that the jets stop in the Galactic halo 
when their radius reaches Rp ~ O.lpc (if nh > 10~* cm~'^) 
and form hot spots. This is illustrated in Fig 1. 

The "decelerating plasmoid model" of Chiang & Der- 
mer (1997) for GRB afterglows makes similar assumptions 
and works out some consequences in more numerical de- 
tail. In particular, in their model the plasmoid expands 
only slowly, and decelerates mainly via sweeping up am- 
bient matter. In their case a plasmoid with a radius of 
about 10"^ to 10"^ pc, ejected by the CGRB with a ki- 
netic energy of 10^" ergs, travels 30 - 100 pc in the first 3 
- 10 days after the GRB (measured in the observer frame) 
and does not slow down to non-relativistic speeds during 
the simulated propagation. 



4. Energy spectrum of CR produced by the plas- 
moids 

4-.1. High-energy limit for the acceleration of CR in the 
plasmoids 

The typical equipartition magnetic fields in the hot 
spots discussed in the previous section may reach B ^ 
(3i?k/-Rp)^^^ ^ 1 G. Synchrotron losses cut off Fermi ac- 
celeration of CR nuclei with mass number A at -E ~ 
TA^Z-^/^{B/G)-^^^ X 10^0 eV. Particle escape cuts off 
Fermi acceleration when the Larmor radius of the ac- 
celerated particles in the plasmoid rest frame becomes 
comparable to the radius of the plasmoid, i.e., above 
E ~ TZiB/G){Rp/0.1pc) X 10^° cV. In the hot spots. 



0.1 pc, r ~ 1 and B ^ 1 G. Consequently, CR with 
E > Z X 10^° eV can no longer be isotropised by accelera- 
tion or deflection in the hot spots. Thus, CR with energies 
above 10^" eV are heavy nuclei in our scenario. 



4-. 2. The energy spectrum at lower energies 

Fermi acceleration in or by the highly relativistic jets 
from GRBs could produce a broken power-law spec- 
trum, dn/dE ^ E^P , with j3 ^ 2.2 below a knee 
around -Eknoo ^ A PeV and (i ^ 2.5 above this en- 
ergy (Dar 1998b). Spectral indices /? ^ 2-3 are expected 
in relativistic shock acceleration with tangled magnetic 
fields (Ballard & Heavens 1992). Galactic magnetic con- 
finement increases the density of Galactic CR by the 
"concentration factor" c — ct]^/Rq, where T]^{E) is the 
mean residence time in the halo of Galactic CR with en- 
ergy E and Rq is the radius of the Galactic magnetic- 
confinement region. With the standard choice for the 
energy dependence of the diffusion constant (observed, 
e.g., in solar-system plasmas (Berezinskii et al. 1990)) one 
gets: T\i ex (_E/Z)^°'^. Consequently, the energy spectrum 
of CR is expected to be 



dn/dE ^ CiE/Ei^neeY 



(4) 



with a ~ (3 + 0.5 ~ 2.7 (~ 3) below (above) the knee. 
The normalisation constant cannot yet be determined - 
except by the rough total-energy argument mentioned in 
Sect. 2 after Eq.(2) - but is universal for all energies. In 
this respect our scenario and the conventional framework 
are in the same status. 

This power-law is expected to continue as long as the 
Galactic magnetic field confines the CRs, i.e., up to ener- 
gies near the "ankle" . 

4-3. Explanation of the "ankle" 

Part of the kinetic energy released by GGRBs is trans- 
ported into the Galactic halo by the jets. Assuming 
equipartition of this energy, without large losses, between 
CR, gas and magnetic fields in the halo during the res- 
idence time of CR there, the magnetic field strength 
Bh in the halo is expected to be comparable to that 
of the disk By, ~ {2LMw[CR]Th/Rl)^/^ ~ 3 ^^G where 
Th ~ 5 X 10^ years is the mean residence time of the bulk 
of the CR in the Galactic halo (see Sect. 4.4) cosmic rays 
with Larmor radius larger than the coherence length A of 
the halo magnetic fields (Lampard et al. 1997), i.e., with 
energy above 



E, 



'ankle 



' 3 X 10^^{ZBij3^G){X/kpc) eV, 



(5) 



escape Galactic trapping. Our "preferred" values indi- 
cated in Eq.(5) are well within the range of possible val- 
ues preferred by other researchers (Lampard et al. 1997). 
Thus, the CR ankle is explained as the energy where 
the mean residence time Th(E) of CR becomes com- 
parable to the free escape time from the halo Tfroo ~ 
1.6(i?h/50 kpc) X 10^ years. This explanation for the an- 
kle was already proposed by one of us (Plaga 1998) in 
the context of a proposed extragalactic origin of hadronic 
CRs, to which the present proposal bears of course some 
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reaches 3 x 10'* at GcV energies (see Fig. 2), correspond- 
ing to a confinement time in the Galactic halo of about 
5 X 10^ years. We interpret the experimentally measured 
"lifetime" of low-energy CR, of the order 2 x 10^ years 
(Berezinskii et al. 1990; Connel et al. 1998), as the resi- 
dence time in the Galactic disk, in which we are located 
as observers, rather than the time since acceleration, as 
in most other models. The mean residence time Td of low- 
energy particles in the disk is given as: 



Td~71i-(dd/dh)'(Al/i^d). 



io'%o"io"io"io"io"io"io"io"io"io'"io" 

E [eV/nucleus] 



Fig. 2. The observed flux of cosmic rays (thin line) as a func- 
tion of primary energy E is well described by a power law that 
changes its slope sharply at only two energies, the "knee" and 
the "ankle". At energies below the ankle it is enhanced (by a 
factor (E/Eanfcie) "■^) over the GGRB source spectrum (thick 
line, a power law with differential power law index of -2.2 be- 
low the knee and ~ -2.5 above it) by way of trapping in the 
Galactic halo magnetic fields 



similarity (in both cases all hadronic CRs have a common 
origin outside the Galactic disk). 

The spectrum of CRs with energies above the ankle that 
do not suffer Galactic magnetic trapping is then the CR 
spectrum produced by the jet, i.e.. 



-2.5. 



dn/dE ^ C{E ankle /Ei^ncc) {E/Eiink\c) , 

E > -Eanklc- 



(6) 



This predicted spectral index of 2.5 has to be compared 
with the one recently determined by the AGASA collab- 
oration (Takeda et al. 1998) of 2.78lo;33. In the conven- 
tional framework - where only UHECRs have an extra- 
galactic origin - the position of the ankle at just the en- 
ergy where Galactic magnetic fields can no longer confine 
CRs is pure coincidence. 

4-. 4- The "knee" in the CR spectrum and low- energy prop- 
agation 

The knee in the CR source spectrum still lacks a com- 
plete explanation (see however Dar (1998b)), but we think 
that a single source population over the whole CR en- 
ergy range offers a more natural framework to understand 
the CR spectrum than the prevailing one where a "fine 
tuning" in the intensity of different source components 
is often required at the knee. The "concentration factor" 
c (see Sect. 4.2 above) over the GGRB source spectrum 



(7) 



If the diffusion coefficient _Dd in the disk is smaller than 
that in the halo (Dh) by about a factor 10^ - due to a 
stronger turbulence in the disk induced by SN explosions 
- then Td is of the observed order of magnitude. 

5. A prediction of our scenario: angular clustering 
of UHE CRs 

Small defiections of UHECR by turbulent magnetic fields 
along their arrival trajectory, with a length d, spread their 
arrival times over 



At - 750(d/50 kpc)2(A/kpc) x 
{ZB/3nGf{E/100 EeY)-^ years. 



(8) 



The nearly isotropic release of CR in the halo by GGRBs 
and the spread in their arrival times, during which many 
GGRBs occur, accounts for the nearly isotropic sky dis- 
tribution of the UHECR. Jets which stop in the Galactic 
bulge or the Galactic disk (see Fig. 1) may produce a de- 
tectable enhancement in the flux of UHECR from these 
directions. 

First data from the Akeno Giant Air Shower Array 
(AGASA) (Hayashida et al. 1996) suggested a clustering 
in the arrival directions of the UHECR, these indications 
have recently been reaffirmed with an enlarged database 
(Takeda et al. 1999). The directions of these clusters do 
not coincide with any known nearby extragalactical ob- 
jects which would seem a plausible candidate for acceler- 
ation the UHECR; this is consistent with a Galactic halo 
origin. Such clustering is predicted if the UHE CRs are 
produced/released in hot spots in the Galactic halo. Their 
accumulated r.m.s. deviation angle over a distance d by 
random walk in the turbulent magnetic field is 



5.4"(d/50 kpc)^/^{\/kpc)^/^{ZB/3^iG) x 

{E/100EeY)-\ 



(9) 



Therefore, if the origin of UHECR is ^ 0.1 pc-sized hot 
spots in the halo, the arrival direction of CR with ener- 
gies around 100 EeV is predicted to cluster around ~ 30 
(arrival time spread x GGRB rate) fixed hot spot posi- 
tions in the sky whose number and spread around them 
decrease with increasing energy. 
An origin of UHE CRs from the decay or annihilation 
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of ultra-massive new particles in the Galactic halo has re- 
cently been discussed (Hillas 1998; Berezinskii et al. 1997; 
Kuz'niin & Tkachev 1998). In such a scenario clustering 
is not possible if the the UHE CRs are nucleons from 
the decay or annihilation of relic particles with mass 
smaller than the Planck mass, mp — y/hc/G w 1.22 x 
10^^ GeV/c , which are distributed isotropically in the 
sky. Even if all the rest mass energy was converted into two 
opposite jets of 10^^ eV nucleons with transverse momen- 
tum of ^ 1 GeV/c with respect to the jet axis, the mean 
number of UHE CR events that such a jet from a typi- 
cal halo distance of 30 kpc could produce in the AGASA 
detector (100 km^) is still only -- 5 x 10"^. 

6. Some experimental data on CRs compared to 
expectations in our scenario 

The CR spectrum which is suggested in our scenario is 
displayed in Fig. 2. As discussed in the previous sec- 
tions it agrees well with the measured spectrum. There 
is no GZK cutoff because of the small distance to the 
sources (0(50 kpc)), over which absorption of UHE nu- 
cleons in the microwave-background radiation is negligi- 
ble. As the major part of the hadronic CRs are not pro- 
duced in SNRs in our scenario, the predicted flux of high- 
energy 7-rays due to proton interaction from these ob- 
jects can be low, thus resolving the second anomaly men- 
tioned in the introduction. A quantitative calculation of 
the expected density of hadronic CRs as a function of 
galactocentric radius in our scenario is difHcult because 
of the ill-understood propagation properties of the plas- 
moids in different part of the Galaxy (halo, bulge, etc.). 
In first approximation this density is predicted to fall with 
a halo scale (i.e., tens of kpc). This is in better agree- 
ment with observations than the ~ 5 kpc scale height of 
SNRs, which is expected to be close to the characteristic 
decay scale in the conventional scenario (third anomaly 
in the Introduction; Strong & Moskalcnko 1998^). Recent 
determinations of the gradient from EGRET data yield 
for the exponential decay scale between 5 kpc and 20 
kpc (read off diagrams in these publications): 16 kpc 
(Hunter et al. 1997), 34 kpc (Erlykin et al. 1996), 23 kpc 
(Strong & Mattox 1996). 

Moreover - due to the near isotropy of their sources - the 
large-scale anisotropy of CRs must be very small in our 
scenario at all energies. This resolves the fourth anomaly 
mentioned in the Introduction. 

7. Properties of decaying plasmoids in the halo 

Electrons suffer large energy losses by synchrotron emis- 
sion and inverse Compton scattering while diffusing out 
of their jet acceleration sites. Thus, CR electrons in the 
Galactic disk could have their origin mainly in Galactic 



^ The propagation models assuming a SNR origin of CRs in 
this publication predict decay scales between 4.5 to 9 kpc. 



SNRs (Ginzburg 1957); this theory has recently received 
experimental support from X-ray (Allen et al. 1997) and 
TeV 7-ray (Tanimori et al. 1998) observations of Cas-A 
and the remnant of SN 1006, respectively. Totani (1998) 
argued that synchrotron radiation from accelerated pro- 
tons is an important emission mechanism in GRB after- 
glows. Assuming hot spots accelerate mainly baryons on 
a time scale of about 1000 years over which they irradiate 
2q50 gj.gg jj^ proton synchrotron radiation, their resulting 
radio intensity will be less than a mJy (Plaga et al. 1999). 
The hot spots are thus not expected to be particularly con- 
spicuous radio sources. The detailed observational conse- 
quences of hot-spot remnants in the halo of the Milky Way 
will be discussed elsewhere (Plaga et al. 1999). 

8. Conclusion 

We proposed an alternative source for the origin of 
hadronic CRs at all energies which is in better accord 
with some observational facts than the conventional one. 
Some of the assumptions on which it is based are certainly 
speculative. However, our ideas seem conservative when 
compared to the audacity of some proposals by outstand- 
ing scholars to explain the puzzling properties of ultra- 
high energy CRs - like, e.g., a violation of Lorentz invari- 
ance (Coleman & Glashow 1998), the existence of a new 
class of super- massive elementary particles in the Galactic 
halo (Berezinskii et al. 1997; Kuz'min & Tkachev 1998; 
Hillas 1998) or supersymmetric light strongly- interacting 
particles (Biermann & Farrar 1998). Further observations 
of GRBs and their afterglows will reveal if the contention 
that a GRB is a common phenomenon in the life of a neu- 
tron star - crucial for our scenario - is tenable. The next 
generation of UHE CR detectors (Mantsch 1996) will rule 
out or confirm our prediction about their angular cluster- 
ing properties. Basic issues of our scenario which require 
more work are the exact nature of the event giving rise 
to a GRB and the observational properties of "hot-spot 
remnants" in the Galactic halo. 
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